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Abstract 

We study the indirect detection of neutralino dark matter using positrons and gamma rays 
from its annihilation in the galactic halo. Considering the HESS data as the spectrum con- 
stituting the gamma-ray background, we compare the prospects for the experiments GLAST 
and PAMELA in a general supergravity framework with non-universal scalar and gaugino 
masses. We show that with a boost factor of about 10, PAMELA will be competitive with 
GLAST for typical NFW cuspy profiles. 
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1 Introduction 



The existence of dark matter in the Universe is generally accepted by the scientific commu- 
nity, though its nature is still unknown. One of the most popular candidates, the lightest 
neutralino, appears naturally in supersymmetric (SUSY) extensions of the standard model. 
It has long been thought that dark matter particles could be observed indirectly by detect- 
ing the products of their annihilations. Such products include gamma-rays, neutrinos and 
anti-matter. Actually, many experiments have been performed or are planned in order to 
detect indirectly (or directly) the presence of such particles in the galactic halo [I]. For ex- 
ample, high-energy positron excess in the cosmic rays provides an opportunity to search for 
the dark matter signal. These fluxes will be measured with unprecedented accuracy by the 
present space-based experiment PAMELA (Payload for Antimatter-Matter Exploration and 
Light-nuclei Astrophysics [2]), and also by future experiments like AMS-02 (Alpha Magnetic 
Spectrometer |3J), designed to be deployed in the International Space Station around 2009-10. 
The precision measurements of the cosmic positron spectrum to be provided by PAMELA will 
be very important to identify signatures of dark matter in our Galaxy. 

On the other hand, the gamma-ray data, richer than the cosmic particle data, also give us 
exciting expectations. The gamma-ray detection experiments have the advantage of being able 
to search for point sources of annihilation radiation. In particular, the center of our Galaxy, 
given its dark matter density, is one of the most promising sources of diffuse gamma-rays 
from dark matter annihilation. In addition to the controversial EGRET data [4], atmospheric 
Cerenkov telescopes like HESS [5j and MAGIC [6j have observed recently a very bright source 
in this direction. Actually, HESS was able to measure in detail the gamma-ray spectrum. 
Although fitting the data with a reasonable SUSY model seems quite difficult, one can use 
them as the astrophysical background for gamma-ray detection [7]. Scheduled to begin its 
five years mission in 2007, GLAST (gamma-ray Large Area Space Telescope [Hj) will be the 
most sensitive gamma-ray telescope in the energy range of interest (1 — 300 GeV). 

Both space-based experiments, PAMELA and GLAST, will provide us their first results 
around 2009 (before the LHC data). The aim of this work is to compare the detection 
prospects of these two experiments in the framework of supergravity (SUGRA) models, taking 
into account the different astrophysical uncertainties: boost factor and halo density profile. 
In particular, the effect of clumps in the galactic halo may increase substantially the positron 
fluxes, and this is parameterised by the so-called boost factor. Likewise, cuspy halo profiles 
may also give rise to larger gamma-ray fluxes. Let us finally remark that although antiproton 
fluxes may be in principle as competitive as the positron ones, we prefer to carry out the 
analysis of antiprotons in detail in a future work |9J. 

The paper is organized then as follows. In section 2 we will review the neutralino charac- 
teristics as a dark matter candidate in SUGRA models. In section 3 we will recall the main 
features of the positron propagation and detection in the light of the PAMELA experiment. 
Section 4 will be devoted to gamma-rays detection and GLAST prospects. In Section 5 we 
will compare in detail these two detection modes. Finally, the conclusions are left for section 
6. 

2 Neutralino dark matter 

For the computation of the positron and gamma-ray fluxes it is crucial to evaluate the 
neutralino annihilation cross section. Of course, for determining this value the theoretical 
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Figure 1: Dominant neutralino annihilation diagrams. Relevant parts of the amplitudes are shown explicitly Terms 
between parenthesis correspond to f u and Z final states in second and fourth diagrams. V and Z are chargino and 
neutralino mixing matrices. 

framework must be established. In particular, we will work in the context of the minimal 
supersymmetric standard model (MSSM). Let us recall that in this framework there are four 
neutralinos, Xi (i = 1)2,3,4), since they are the physical superpositions of the fermionic 
partners of the neutral electroweak gauge bosons, called bino (-B ) and wino (W3), and of 
the fermionic partners of the neutral Higgs bosons, called Higgsinos (H®, H^). Thus one can 
express the lightest neutralino as 

X? = Z n B° + Z 12 Wl + Z n H° d + Z lA H Q u . (2.1) 

It is commonly defined that Xi is mostly gaugino-like if P = \Zu\ 2 + \Z\ 2 \ 2 > 0.9, Higgsino-like 
if P < 0.1, and mixed otherwise. 

In figure[I]we show the relevant Feynman diagrams contributing to neutralino annihilation. 
As was remarked e.g. in Refs. [TQl CD], the annihilation cross section can be significantly 
enhanced depending on the SUSY model under consideration. We will concentrate here on 
the SUGRA scenario, where the soft terms are determined at the unification scale, Mqut ~ 
2 x 10 16 GeV, after SUSY breaking, and radiative electroweak symmetry breaking is imposed. 

2.1 Supergravity models 

Let us discuss first the minimal supergravity (mSUGRA) scenario, where the soft terms of 
the MSSM are assumed to be universal at Mqut- Recall that in mSUGRA one has only four 
free parameters: the soft scalar mass m, the soft gaugino mass M, the soft trilinear coupling 
A, and the ratio of the Higgs vacuum expectation values, tan/3 = (H^}/(H d ). In addition, 
the sign of the Higgsino mass parameter, fx, remains also undetermined by the minimization 
of the Higgs potential. 

Since in this scenario the lightest neutralino Xi is mainly bino, only Z\\ is large and then 
the contribution of diagrams in Fig. [1] will be generically small, the first being suppressed by 
/ masses. As a consequence, for example to reproduce in mSUGRA the present experimental 
accesible regions (EGRET data) is not possible pj2] for a typical halo model. 

However, as discussed in detail in Ref. (TOl CE2] in the context of indirect detection, the 
annihilation cross section can be increased in different ways when the structure of mSUGRA 
for the soft terms is abandoned. In particular, it is possible to enhance the annihilation 
channels involving the exchange of the CP-odd Higgs, A, by reducing the Higgs mass. In 
addition, it is also possible to increase the Higgsino components of the lightest neutralino 
■^13,14- Thus annihilation channels through Higgs exchange become more important than in 
mSUGRA. This is also the case for Z, xf, an d Xi-exchange channels. As a consequence, 
positron and gamma-ray fluxes will be increased. 
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In particular, the most important effects are produced by the non-universality of Higgs 
and gaugino masses. These can be parameterised, at Mqut, as follows 

m 2 Hd = m\l + 5 1 ) , m 2 Hu =m 2 {\ + 5 2 ) , (2.2) 

and 

M X =M , M 2 = M (1 + 4) , M 3 = M(l + 6' 3 ) . (2.3) 
We will concentrate in our analysis on the following representative cases: 
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Case a) corresponds to mSUGRA with universal soft terms, cases b), c) and d) correspond 
to non-universal Higgs masses, and finally cases e) and f) to non- universal gaugino masses. 
The cases b), c), d), and e) were discussed in Ref. [10 j . and are known to produce gamma-ray 
fluxes larger than in mSUGRA, whereas case f) enhances efficiently annihilation processes 
because of the pure wino nature of the LSP. 

Let us remark that for the evaluation of the gamma-ray fluxes and the positron fluxes 
at production we use the last DarkSusy released version [13J. Concerning the computation 
of the fluxes observed in the solar neighborhood in the case of adiabatically compressed ha- 
los we use our own code. To solve the renormalization group equations (RGEs) for the soft 
SUSY-breaking terms between Mqut and the electroweak scale, we use the Fortran package 
SUSPECT [14]. Concerning the positron spectrum observed by PAMELA, we use the anni- 
hilation cross section given by DarkSusy at the source, but solve ourselves the propagation 
equation and calculation of the flux measured on the Earth. 



2.2 Experimental and astrophysical constraints 

We have taken into account in the computation several experimental and astrophysical bounds. 
These may produce important constraints on the parameter space of SUGRA models, restrict- 
ing therefore the regions where the gamma-ray and positron fluxes have to be analyzed. 

In particular, concerning the astrophysical constraints, the bounds 0.1 < Odm^ 2 ^ 0.3, 
on the relic neutralino density computation has been considered. Due to its relevance, the 
WMAP narrow range 0.094 < Q.uMh 2 < 0.129 has also been analyzed in detail. 

Concerning the experimental constraints, the lower bounds on the masses of SUSY particles 
and on the lightest Higgs have been implemented, as well as the experimental bounds on the 
branching ratio of the b — > process and on a^ USY . Note that we are using \i > 0. We will 
not consider in the calculation the opposite sign of \x because this would produce a negative 
contribution for the g^ — 2, and, as will be discussed below, we are mainly interested in positive 
contributions. Recall that the sign of the contribution is basically given by fJ,M 2 , and that 
M, and therefore M 2 , can always be made positive after performing an U(1)r rotation. For 
a^ USY , we have taken into account the recent experimental result for the muon anomalous 
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magnetic moment [E], as well as the most recent theoretical evaluations of the standard 
model contributions [IB] . It is found that when e + e~ data are used the experimental excess in 
(g^ - 2) would constrain a possible SUSY contribution to be a^ USY = (27.1 ± 10) x 1CT 10 . At 
2a level this implies 7.1 x 1CP 10 < a^ USY < 47.1 x 10~ 10 . It is worth noticing here that when 
tau data are used a smaller discrepancy with the experimental measurement is found. In order 
not to exclude the latter possibility we will discuss the relevant value a^ USY = 7.1 x 1CT 10 . 

On the other hand, the measurements of B — > X s ^ decays at CLEO p2] and BELLE [18J, 
lead to bounds on the branching ratio b — * sj, In particular we impose on our computation 
2.33 x 10~ 4 < BR(b — > 57) < 4.15 x 10~ 4 , where the evaluation is carried out using the routine 
provided by the program micrOMEGAs [19]. This program is also used for our evaluation of 
a^ USY and relic neutralino density. 



3 The positrons in the Galaxy 
3.1 The background 

The conventional sources of cosmic rays are believed to be supernovae and supernovae rem- 
nants, pulsars, compact objects in close binary systems, and stellar winds. Observations of 
synchrotron emission and 7— rays reveal the presence of energetic particles in this objects, thus 
testifying to efficient acceleration processes. Propagation in the interstellar medium changes 
the initial composition and spectra of cosmic ray species due to the energy losses (ionization, 
Coulomb scattering, bremsstrahlung, inverse Compton scattering, and synchrotron emission), 
energy gain (diffusive re-acceleration), and other processes (i.e., diffusion and convection by 
the galactic wind) [20J. The destruction of primary nuclei via spallation gives rise to secondary 
nuclei and isotopes which are rare in nature (i.e., Li, Be, B), antiprotons, and pions (ir^, ir°) 
that decay producing secondary and 7— rays. The abundance of stable (Li, Be, B, Sc, Ti, 
V) and radioactive (Be 10 , Al 2G , CI 36 , Mn 54 ) secondaries in cosmic rays are used to derive 
the diffusion coefficient and the halo size [21] . The result of recent simulations agrees with 
measurements of the low-energy positron flux in the cosmic rays [21]. The fitting functions 
for high energy positrons, primary electrons, and secondary electrons have been constructed 



5, ?'° (£) = 1 + i^S^ (GeV- W VV-') , (3.5) 

< £ > - i + im^Tmm WV'sr-) , (3.6) 

C e + C (£) = ^ (GeV- 1 cm- 2 s- 1 sr- 1 ) . (3.7) 

e y ' 1 + 650£ 2 - 3 + 1500£ 4 - 2 v ' v ' 

3.2 The propagation of a positron 

Concerning the propagation of positrons, we will use the "diffusion model" in which the 
random walk is described by the diffusion equation, 

^f(E : r)=K(E)Af(E : r) + -^[b(E)f(E,r)]+Q(E,r) , (3.8) 



6 



where f(E, r) denotes the number density of particles per unit of volume and energy and 
Q(E,r) is the source (positron injection) term generated by neutralino annihilation. The flux 
of positrons with high energy (E 3> m e ) around the Sun is given by: 

^f(E,r Q )- (3.9) 

Above a few GeV, positron energy losses are dominated by synchrotron radiation in the galactic 
fields and by inverse Compton scattering on stellar light and on CMB photons. The energy 
loss rate b(E) depends on the positron energy E through : 

HE) = £- , (3.10) 

where we have set the energy reference Eq to 1 GeV and the typical energy loss time is 
te = 10 16 s. We have also assume that the space diffusion coefficient K(E) is written on the 
form 

K{E) = kJ^X , (3.11) 

where the diffusion coefficient at 1 GeV is Kq = 3 x 10 27 cm 2 s _1 with a spectral index of 
a = 0.6 [22]. 

To solve Eq. (|3.8p we make the hypothesis that the positron fluxes coming from dark 



matter annihilation are in equilibrium in the present universe i.e., ^f(E,r) = 0. Thus 
Eq. (13. 8p simplifies into 



K e a Af(e,r) + ^- 



—f(e,r) 

TE 



+ Q(e,r) = 0, (3.12) 



where we have defined e = E/Eq. 

Concerning the boundary condition, we impose the "free escape" one. This means that the 
positron density drops to zero on the surface of the diffusion zone. Moreover, we neglect the 
positrons coming from outside our Galaxy. The diffusion zone (Galaxy) will be parameterised 
by a cylinder of half height L = 3 kpc and radius R = 20 kpc. We will see afterwards that 
this boundary condition does not affect too much the result around the Sun. 

We can see from Eq. (|3.8p that the high-energy positrons mainly come from a region 
within few kpc from the solar system. In fact, this depends mainly on the parameters K 
and b appearing there. Indeed, positrons far from the Earth lose their energies during the 
propagation, and consequently they contribute to the low-energy part of the spectrum. Using 
Eq. (|3.8p we can easily approximate the distance in which positrons travel without significant 
energy loss: 



r « y/ = 1.8 x (E/lGeV)- - 2 kpc , (3.13) 

which gives r ~ 1.8 kpc for 1 GeV and r ~ 0.72 kpc for a 100 GeV positron. This implies 
that the influence of the different kinds of dark matter density profiles on the positron fluxes 
is negligible (recall that all simulations give rise to the same kind of profile around the solar 
system). Thus the astrophysical dependence of the indirect detection of dark matter through 
positrons, instead of being similar to the case of the indirect detection through gamma-rays, 
looks more like the case of the direct detection of dark matter. 
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E e+ (GeV) 

Figure 2: The positron flux from neutralino annihilation in the non-universal case /) M 2 = 0.5M with tan/5 = 35, 
fj, = 0, m = 1 TeV, M = 1 TeV (corresponding to a neutralino mass m x — 401) as a function of positron energy E 
before (dashed-line) and after (continuous line) propagation in the interstellar medium. 



To solve Eq. (13. 12ft we use the treatment given in Ref . [23] , where the solution is expressed 
with the Green function G as 0: 



f(e,r) 



de s / d 6 r s G(e,e s ,r,r s ) Q(e s ,r s 



(3.14) 



with 

G(e,e s ,r, r a ) 
and 
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r=[x 2 + y 2 + z 2 ] 1/2 , p =l( x - Xs f + (y-y s fyi\ (3.16) 

The expression for V(e, e s , z, z s ) is obtained using the original method of "electrical images' 
firstly introduced by Blatz and Edsjo [22] 



x V(e,e s ,z,z s ) , (3.15) 
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<t>'n( z s)<t>n( z )} 



(3.17) 



where X n = Kq x (n — 1/2)^, = Kq x n^, and n (z) = sin£; n (L— \z\), <fi' n (z) = 
sink' n (L — z), with k n = 1/L(n — l/2)ir, k' n = 1/Lnn. We have checked that (as was al- 
ready noticed in [23]) usually just a few eigenfunctions 4> n , <f/ n , need to be considered for the 
sum in Eq. (I3.17P to converge. We illustrate the effect of the propagation on the spectrum in 
Fig. [2l Clearly, this is softened after the propagation through the interstellar medium. 

3.3 The boost factor 

Recently, many works based on N-body simulations discussed the effect of large inhomogeneity 
(clumps) in the galactic halo. In particular, it has been shown that dark matter annihilation 



!An interesting alternative to solve the propagation equation is given in [24J. 
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tanp= 35 ; (m, M)= (1 TeV, 1 TeV) ;^ = 0.5M 




E (GeV) 

e+ 

Figure 3: Positron flux dependence on the dark matter halo profile. A non-cuspy isothermal profile (continuous 
line) and the very cuspy Moore et al. profile with adiabatic compression (dashed line) are shown. Black upper 
curves show the background, and the signal plus background. 



and detection rates may be largely enhanced around the solar system [25j [26], [23] • The effect 
is parameterised by a multiplicative factor, the so-called boost factor (BF), which is defined 
by the ratio of the signal fluxes with inhomogeneity and without inhomogeneity, 

fy d 3 xp clu f v d i xp dum 
BF= 2 ^~ 2T7 — el , 3.18 

JV d x Psmooth P0 V 

where the region of integration is taken to be V ~ (a few kpc) 3 around the solar system. 
Obviously, the boost factor is one only when the density p is constant, otherwise it is always 
larger than one. Clustering scenarios suggest a boost factor of ~ 2 — qj. 



3.4 Profile dependence 

We show in Fig. [3]the influence of the dark matter distribution on the positron fluxes measured 
on the Earth. We have calculated the positron flux in two extreme distributions around the 
galactic center: a non-divergent isothermal profile and a cuspy one (Moore et al. profile with 
adiabatic compression [12]). Our numerical result clearly agrees with the approximated result 
discussed in Eq. (|3.13h . confirming that most of the positrons observed would be produced 
in the vicinity of the solar system. We remark also that a Moore et al. compressed profile 
affects mainly the low-energy positrons in the spectrum. Indeed, a more cuspy profile will 
enhance the flux of positrons coming from the galactic center, i.e. far away from the solar 
neighborhood. They will lose considerable amount of energy during their propagation. As a 
consequence, such annihilation process will contribute to the low-energy part of the positron 
spectrum, as we can see in Fig. [3l 

2 Let us remark however, that in a recent study a model where the boost factor depends on the energy has been 
proposed [23]. Although a detailed analysis of this possibility is lengthy, discussions with one of the authors of 
Ref. [23] lead us to believe that, due to the characteristics of our positron spectrum, the results of our analysis 
would not be crucially affected. 
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3.5 The PAMELA experiment 



The PAMELA experiment has been launched recently into space. One of the PAMELA 
primary goals is the measurement of the cosmic positron spectrum up to 270 GeV [2]. The 
geometric acceptance in the standard trigger configuration is ~ 20 cm 2 sr. We consider in our 
analysis a three years mission and assume Gaussian statistics to determine the \ 2 ■ 

Let us explain now how we can define that a positron signal arriving from the galactic halo 
will be distinguished effectively as a "signal" from the point of view of an experiment. We will 
mainly follow in this study the work by Lionetto, Morselli and Zdravkovic [27] and Hooper 
and Silk [28]. Let us call (j) susy the signal, 4> bkg the background, and 4> tot = cb susy + (j) bk9 the 
total flux which will be observed by PAMELA. We will divide our energy range between 1 
to 500 GeV in 20 energy bins (N=20) logarithmically evenly spaced. For the discrimination 
between the signal and the background, we calculate the reduced \ 2 ■, X 2 e d : 

5 _ 1 \^ \E n <P n ~ ^n'Pn ) 
Xred-Jj^ 7^)2 > ^.19) 

n=0 

where a n is the error, assuming Gaussian statistic, on the measured value of the flux multiplied 
by E\. For a "discovery signal", we request the condition^] x 2 e d > 1- 
It is easy to check that this a n , for a given E n , can be expressed as : 




(3.20) 



where A is the geometrical factor or acceptance of the experiment mentioned above, and T 
is the exposure time (3 years in our study). This means that a signal at one a in one energy 
bin E n will give us N hits in the range N ± viV after T seconds (Gaussian law). 

To illustrate the discussion we show in Fig. [4] the signal that we can expect for two 
different points in the non-universal SUGRA parameter space (in particular for the case 
/) Mi = 0.5M in Eq. (12. 4)1 ). Let us remark that in addition to the propagation effects, as 
positrons approach the solar system, their interactions with the solar wind and magnetosphere 
can become important. These effects are called solar modulation. Alternatively, if one assumes 
that the effects of solar modulation are charge-sign independent, their impact can be removed 
considering the ratio of positrons to positrons plus electrons at a given energy rather than the 
positron flux alone. This quantity, the so-called positron fraction is often used to minimize 
the uncertainties associated with modelling the impact of solar modulation. As we can convert 
the fluxes to positron fractions by using fluxes of background positrons and electrons given 
by the Eqs. (|3.5I) - (I3.7I) . we can replace the number of events used in Eq. (13 . 19|) with ratios 
of positrons to positrons+electrons to reduce the effects of solar modulation in our results. 
Because there are substantially more electrons than positrons observed, we can assume that 
there are negligible errors associated with the electron flux. 

3 It is worth noticing here that Xrcd > 1 does not really mean a la discovery in our limit. To be coherent, the 
important point is to compare two experiments with the same criterion of discrimination. Thus Xrcd > 1 wm De 
the criterion of comparison common to PAMELA and GLAST in our study. 
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E + (GeV) E + (GeV) 

Figure 4: PAMELA expectations for positrons in the non-universal case /) M 2 = 0.5M with tan/3 = 35, A = 0, 
m = 1 TeV, M = 1 TeV (left) and m = 2 TeV M = 2 TeV (right), corresponding to a neutralino mass of m x = 401 
GeV and m x — 820 GeV, respectively. The boost factor BF has been set to 5 in an isothermal profile. The error 
bars shown are those projected for the PAMELA experiment after 3 years of observation. The lower (blue) line is 
the background-only prediction. The reduced \ 2 ls 312 in the first case giving a signal clearly distinguishable from 
the background whereas x 2 e d ls on ^y 1-9 m the second case and gives a good fit of the background. 





a (kpc) a (3 7 J(l(T 5 sr) 


NFW 
NFW C 
Moore et al. 
Moore c 


20 1 3 1 1.264 x 10 4 
20 0.8 2.7 1.45 1.205 x 10 7 
28 1.5 3 1.5 5.531 x 10 6 
28 0.8 2.7 1.65 5.262 x 10 8 



Table 1: NFW and Moore et al. density profiles without and with adiabatic compression (NFW C 
and Moore c respectively) with the corresponding parameters, and values of J(AO). 



4 The gamma-rays in the Galaxy 
4.1 The flux 

The spectrum of gamma-rays generated in dark matter annihilation and coming from a di- 
rection forming an angle tp with respect to the galactic center is 

dN" 1 1 r 

where the discrete sum is over all dark matter annihilation channels, dN^/dE^ is the dif- 
ferential gamma-ray yield, (a{v) is the annihilation cross section averaged over its velocity 
distribution, and p is the dark matter density. 



It is customary to rewrite Eq. (|4.21|) introducing the dimensionless quantity J (which 
depends only on the dark matter distribution): 



8.5 kpc V0.3 GeV/cm J 



J W = 771Z: L 3 ^v/^ / p 2 {r{lA))dl. (4.22) 



line of sight 
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After having averaged over a solid angle, Af2, the gamma-ray flux can now be expressed as 



$ 7 (£ 7 ) = 0.94 x 1(T 13 cm~ 2 s" 1 GeV" 1 ™ _1 



sr 



^dN, % ( ( ai v) \ /100 GeV\% /Af , x A „ /4nn . 

The value of J(AO)AO depends crucially on the dark matter distribution. The different 
profiles that have been proposed in the literature can be parameterised as 

p(r) - PQ[l + (W]^ )/a (4 24) 

where po is the local (solar neighborhood) halo density and a is a characteristic length. Al- 
though we will use po = 0.3 GeV/cm 3 throughout the paper, since this is just a scaling factor 
in the analysis, modifications to its value can be straightforwardly taken into account in the 
results. N-body simulations suggest a cuspy inner region of dark matter halo with a distribu- 
tion where 7 generally lies in the range 1 (NFW profile [29]) to 1.5 (Moore et al. profile [30] fl, 
producing a profile with a behaviour p(r) oc r~ 7 at small distances. Over a solid angle of 
10~ 5 sr, such profiles can lead to J (AO) ~ 10 4 to 10 7 . Moreover, if we take into account the 
baryon distribution in the Galaxy, we can predict even more cuspy profiles with 7 in the range 
1.45 to 1.65 (J(Af2) ~ 10 7 — 10 8 ) through the adiabatic compression process (see the study of 
Refs. [32, [12]). We summarize the parameters used in our study and the values of J for each 
profile in table[Q It is worth noticing here that we are neglecting the effect of dumpiness, even 
though other studies showed that, depending upon assumptions on the clumps distribution, 
in principle an enhacement of a factor 2 to 10 is possible [26]. Thus our predictions below for 
GLAST of the gamma-ray flux from the galactic center from dark matter pair annihilations, 
are conservative in this respect. 



4.2 The background 

HESS [5] has recently measured the gamma-ray spectrum from the galactic center in the 
range of energy ~ [160 GeV-10 TeV]. The collaboration claim that the data are fitted by a 
power-law F(E) = F E^ Y , with a spectral index a = 2.21 ± 0.09 and F = (2.50 ± 0.21) x 
10 cm" - 2 s- 1 TeV' 1 The data were taken during the second phase of the measurement 
(July-August, 2003) with a x 2 of 0.6 per degree of freedom. Because of the constant slope 
power-law observed by HESS, it results possible but difficult to conciliate such a spectrum 
with a signal from dark matter annihilation [T2], [33] . Indeed, final particles (quarks, leptons 
or gauge bosons) produced through annihilations give rise to an spectrum with a continuously 
changing slope. Several astrophysical models have been proposed in order to match the HESS 
data [34]. In the present study we consider the astrophysical background for gamma-ray 
detection as the one extrapolated from the HESS data with a continuous power-law over 
the energy range of interest (« 1 - 300 GeV). As was recently underlined in 0, GLAST 
sensitivity is affected by the presence of such an astrophysical source. Note that neutralino 
masses obtained in our parameter space < 1 TeV avoid any conflict with the observations of 
HESS. 

In addition, we have also taken into account the EGRET data [4] in our background at 
energies below 10 GeV, as they can affect the sensitivity of the analysis. Indeed, the extrapo- 
lation of the gamma-ray fluxes measured by HESS down to energies as low as 1 GeV is likely 

4 For a different viewpoint favoring cored distributions, see Ref. |31| . 
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Figure 5: GLAST expectations for gamma-ray fluxes in the non-universal case /) M 2 = 0.5M with tan/3 = 35, 
A = 0, m = 500 GeV, M = 500 GeV (left) and m — 2 TeV M = 2 TeV (right), corresponding to a neutralino mass 
of m x — 185.4 GeV and 820 GeV, respectively. A NFW profile is used in both cases. The error bars shown are those 
projected for the GLAST experiment after 3 years of observation. The lower (blue) line is the background-only 
prediction. The reduced x 2 is 134 in the first case giving a signal clearly distinguishable from the background, 
whereas it is only 0.02 in the second case and gives a very good fit of the background. 



to be an underestimation of the gamma-ray background in the galactic center, as EGRET 
measurements are one to two orders of magnitude higher than the HESS extrapolation. We 
have decided then to take as our background an interpolation between the HESS extrapola- 
tion and the EGRET data below 10 GeV to stay as conservative as possible in evaluating the 
gamma-ray background. 

4.3 The GLAST experiment 

The space-based gamma-ray telescope GLAST [8] is scheduled for launch in 2007 for a five 
years mission. It will perform an all-sky survey covering a large energy range 300 GeV). 

With an effective area and angular resolution on the order of 10 4 cm 2 and 0.1° (AS7 ~ 1CT 5 
sr) respectively, GLAST will be able to point and analyze the inner center of the Milky Way 
(~ 7 P c )- 

In Fig. [5] we show the ability of GLAST to identify a signal from dark matter annihilation 
in the non-universal case /) M 2 = 0.5M with tan (3 = 35 for (m = 500 GeV, M = 500 GeV) 
and (m = 2 TeV, M = 2 TeV) giving a neutralino mass of 185.4 GeV and 820 GeV respectively. 
Concerning the requested condition on the xf-ed f° r a signal discovery, we have used an analysis 
similar to the one considered in the case of positron detection with the PAMELA experiment 
in section [331 with a three years mission run. The error bars shown are projected assuming 
Gaussian statistic, and we adopt a power-law background extrapolated from the HESS data. 
We find that with a Xred °f 1^4, GLAST would potentially be able to detect dark matter 
annihilation radiation from a neutralino of mass 185.4 GeV, whereas with a Xred °f 0-02, a 
signal coming from a neutralino of mass 820 GeV will be below its sensitivity. 
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5 Discussion 



We have plotted in Fig. [6] the isocontour x 2 re d = f° r PAMELA and GLAST, each after three 
years of observations, in the six cases of Eq. (12.41) . and in the SUGRA parameter space (m, M) 
for tan/3 = 35, A = and \x > 0. We have calculated it for boost factors between 1 and 10 4 , 
and three typical halo profiles, isothermal (ISO), NFW, and NFW with adiabatic compression 
(NFW C ). We have taken into account in the computation the astrophysical and experimental 
constraints discussed in section 2.2. In particular, the brown-dark region is excluded by the 
experimental constraints and the stau as the lightest SUSY particle. The region to the right 
of the black dashed line corresponds to a^ USY < 7.1 x 10~ 10 , and would be excluded by e + e~ 
data. The region between thick and thin solid contours fulfils 0.1 < — ^ ^ (^ ne thick 

contour indicates the WMAP range 0.094 < Q^oh 2 < 0.129). Notice however that we did not 
rescale the fluxes according to the cosmological abundance, since this procedure would affect 
photon and positron in the same way and our work focuses on the comparison between the 
two kinds of detection. As usually discussed in the literature, for points with relic density 
away from WMAP values, scenarios with non-thermal production or dilution can be invoked 
to accomodate the abundance of neutralinos. 

We clearly see in the figure that the parameter space of mSUGRA will be reachable only 
in astrophysical scenarios where the halo profile is very cuspy, like NFW C (recall that this 
profile is also very similar to a Moore et al. profile without compression), or very clumpy with 
BF > 100. Notice in this sense that for a given BF (halo profile), the region to the right of 
that blue-dashed (red-dotted) line will correspond to Xred < 1> an d therefore no observation. 
These results are in accordance with a similar study made in the framework of mSUGRA 
for the antiproton detection [27]. In case b), non-universality in m# 2 (62 > 0) increase the 
Higgsino fraction of the neutralino: the annihilation process through Z exchange is enhanced, 
permitting observations in NFW-type profiles or BF ~ 10. Cases c) and d) are similar in the 
sense that the non-universality for niH 1 {Si < 0) open a broad region of A— pole, increasing 

the gamma ray and positron fluxes through XX — > ff, favored in high tan/3 regimes. Profiles 
even less cuspy than NFW and boost factor < 10 could be sufficient for an observation by 
GLAST or PAMELA. The phenomenon is even more spectacular in case e) where the non- 
universality in the gluino sector, M3 (5' 3 < 0), acts on /i and tua through the RGEs. As 
a consequence, the A— pole region is broader and allows observations even for BF ~ 1. In 
case f), the neutralino in mainly Wino and then annihilation to positrons is favored (through 
final states). In this scenario, a clumpy profile with a boost factor of about 1 would be 
comparable with a cuspy NFW profile. Unfortunatly, the amount of relic density in either 
cases is not sufficient to account for the WMAP results. As a conclusion, we can say that 
generically, GLAST and PAMELA will have similar discovery prospects for a NFW profile 
and BF ~ 10, or for a NFW C profile and BF ~ 1000. Obviously, a boost factor as large as 
the latter is unrealistic, and in the case of halo models with adiabatic compression PAMELA 
cannot really be competitive with GLAST. 

6 Conclusions 

We have studied the indirect detection of neutralino dark matter using positrons and gamma 
rays from its annihilation in the galactic halo. Considering the HESS data as the spectrum 
constituting the gamma-ray background, we have compared the prospects for the experiments 
GLAST and PAMELA in a general SUGRA framework with non-universal soft gaugino and 
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Figure 6: The isocontour xl ed = 1 f° r PAMELA (blue-dashed lines) and GLAST (red-dotted lines) for different 



scalar masses. We have shown that with a dumpiness (boost) factor of about 10, PAMELA 
will be competitive with GLAST for a typical NFW cuspy profile. However, for the case of a 
NFW compressed profile, PAMELA cannot be competitive with GLAST since an unrealistic 
dumpiness factor of about 1000 would be necessary. For the future, it would be interesting 
to compile all the detection modes (antiproton, deuteron, positron and gamma rays) in order 
to carry out a complete analysis of the parameter space in a general SUGRA framework [9]. 
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